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We introduce with geometric means a density matrix decomposition of a multipartite quantum 
system of a finite dimension into two density matrices: a separable one, also known as the best 
separable approximation, and an essentially entangled one, which contains no product states com¬ 
ponents. We show that this convex decomposition solving the separability problem, can be achieved 
in practice with the help of an algorithm based on linear programming, which in the general case 
scales polynomially with the dimension of the multipartite system. Furthermore, we suggest meth¬ 
ods for analyzing the multipartite entanglement content of the essentially entangled component and 
derive analytically an upper bound for its rank. We illustrate the algorithm at an example of a com¬ 
posed system of total dimension 12 undergoing loss of coherence due to classical noise and we trace 
the time evolution of its essentially entangled component. We suggest a “geometric” description of 
entanglement dynamics and show how it explains the well-known phenomena of sudden death and 
revival of multipartite entanglement. 
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I. INTRODUCTION 

Though quantum entanglement is a concept which has 
attracted much of the attention of physicists working in 
various fields [l|, still, there remains room for further 
progress on its understanding 0. One of the main open 
problems is the efficient detection and characterization 
of multipartite entanglement of density matrices repre¬ 
senting open quantum systems undergoing non-unitary 
evolution Q. 

All experimentally addressable information about a 
quantum physical system is summarized in its density 
matrix p. We focus on a multipartite quantum system, 
which comprises a finite number AT < oo of parts A4 
numerated by index k = 1,...,K, each of which has 
the Hilbert space of a finite dimensionality Nk, whence 
Y\k=i = N is the dimensionality of the Hilbert space 
of the entire system. This system -assembly of parts, is 
called entangled (or inseparable) if and only if its density 
matrix cannot be caste as a statistical sum 
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of pure states |q!^) of each part. This condition provides 
the most general case of entangled systems opposite to 
a separable quantum system comprised of statistically 
independent element^ where Eq. o holds as an equality. 

Many approaches 0 have been developed so far aiming 
to answer the question whether or not a density matrix 


is separable. Concerning exact analytic results, up to 
now, there is no method applicable to the multipartite 
problem, and we believe that such a solution does not 
exist at all. An algorithmic solution to the “decision” 
problem Q associated with separability has been con¬ 
jectured to be a NP hard problem but valuable progress 
has been done (mainly on the bi-separability problem) 
in approaches 0-0 where semidefinite programming is 
merged with analytic criteria m- 

In this work we provide a geometric point of view on 
the problem of inseparability that suggests an efficient 
solution based on linear programming. Employing simple 
geometric arguments we suggest an algorithm that results 
to a unique decomposition of the density matrix as 

P = (1 - B)pse-p + Bpent (2) 

where psep is, what we call in this work, the separable 
component, pent the essentially entangled part which can¬ 
not have any separable states as components and H is a 
positive number in the range [0,1]. Obviously, the de¬ 
composition, Eq. implies that the state p is separable 
in all K parts only for H = 0. 

The decomposition in Eq.([2]) has been initially intro¬ 
duced in [l^l without resorting to a geometric picture and 
the component (1 — B)psep is widely known as the best 
separable approximation of the density matrix p. In that 
same seminal work, the uniqueness of the decomposition 
has been proven for the multipartite case and a strict 
upper bound on the rank of the component pent for the 
biseparable case. In this work we generalize the latter 
to the multipartite case, proving that the rank of pent is 
upper bounded by a number related to the dimensions of 
the total system and those of the sub-elements. 


2 


On a practical level, we show that linear programming 
algorithm combined with a simple optimization tech¬ 
nique allows one to efficiently find the decomposition of 
a generic density matrix 
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When this limit is reached, the decomposition in Eq.Q 
yields Eq.© with B = (J2Zi 
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entangled states 

Is known that the linear programming method, in the 
general case, scales polynomially with the dimension of 
the vector space where it is applied. Employing the fact 
that M + m < N'^ in Eq.©, where N is dimension of 
quantum assembly under study, we show that the pro¬ 
posed algorithm yielding the decomposition Eq.© scales 
as (2iV4)3. 

The paper is structured as follows. In Section El we 
introduce the idea of the decomposition Eq.© and il¬ 
lustrate its uniqueness with a simple geometric picture 
generalizing the Bloch vector representation of a two- 
level system. This picture also helps us to analyze some 
properties of pent and we conclude this section with a the¬ 
orem setting an upper limit on its rank. In Section Hill we 
present a version of an efficient linear programming algo¬ 
rithm allowing one to explicitly find the decomposition 
Eqs.©-© for a generic density matrix. In Section ITVl 
we suggest ideas for characterizing entanglement of the 
component pent which naturally reflects the entanglement 
properties of p. In Section |V] we present a physical exam¬ 
ple which demonstrates the implementation of the tech¬ 
nique introduced in previous sections and connects them 
with known notions in open quantum system dynamics. 
We conclude by the discussion in Section |VT1 


II. THE GEOMETRIC IDEA OF 
DECOMPOSITION AND PROPERTIES OF THE 
ESSENTIALLY ENTANGLED PART 

All possible density matrices of a quantum system with 
a Hilbert space of dimension iV, are comprising a convex 
set of positive Hermitian matrices of unit trace. This 
set can be viewed as a manifold in the vector space of 
Hermitian matrices endowed with the a metric - the 
Hilbert-Schmidt inner product tr[pi,pj]. The require¬ 
ment of the unit trace in this representation means that 
the inner product of a vector representing a density ma¬ 
trix and a vector representing the unit matrix equals to 
unity. Henceforth we call this manifold “Liouville vector 
space”. Furthermore, the density matrix of a pure state 
has rank one, which implies that the length of the vec¬ 
tor corresponding to a pure state, equals to unity. The 
density matrix manifold is thus a convex hull at the unit- 
length length vectors having unit projection on the unity 
matrix. 


A natural basis exists for such a vector space suggested 
by the properly normalized generators of the uni¬ 
tary group SU{N), including the unity I = VNgQ . This 
basis allows one to cast a NxN density matrix of a quan¬ 
tum system as p = with = Trlgf ^the 

iV^ real vector components. This geometric picture is in 
direct analogy to the Bloch vector for two-level systems. 
The pure quantum states lay at the surface of the unit 
hypersphere, Tr[jO^] = contrast to 

the Bloch vector of 2-dimensional pure quantum states, 
these states do not cover all the surface of the hyper¬ 
sphere of dimension N'^ — 2 but are confined at a manifold 
of lower dimensionality, N'^ — N. This can be easily un¬ 
derstood when the characteristic polynomial Det[A —^ = 
+ ci({rJ)A^“^ -f C 2 ({ri})A ^“2 -f ... -k CAr({rJ) of a 
pure state is considered. The unit trace requirement en¬ 
sures that Ci{{ri}) = —1, while the rank 1 requirement 
implies the constraints Cm{{ri}) = 0 for m = 2,... ,N. 

The set of N conditions on the N'^ components of 
the vector representing a pure state, constrains the vec¬ 
tor to lay on a restricted manifold of lower dimension 
(A^^ — a) at the surface of the unit hypersphere. As 
a consequence, the density matrices for quantum sys¬ 
tems of dimension A > 2 do not ‘fill’ the whole inner 
part of the unit hypersphere, but they are laying inside 
an (A^ — 2)-dimensional body formed as a convex hull 
of the pure states of the (A^ — A) dimensional mani¬ 
fold. This convex hull plays the role of the Bloch ball 
for higher dimensions of the Hilbert space and has been 
exhaustively studied in [l^ for the case of 3-dimensional 
systems. The convex hull is touching the unit hyper¬ 
sphere only for the pure states while its outer hypersur¬ 
face, which we denote by Sch-: is naturally the border 
between positive and non-positive Hermitian matrices of 
unit trace. Therefore Sch consists only of the degenerate 
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FIG. 1: A symbolic illustration of the geometric structure of density matrices and of the decomposition Eq. ©. 


density matrices which have at least one zero eigenvalue. 
In Fig[T](b) we symbolically illustrate the convex hull of 
pure states, such that all density matrices are inside this 
body. 

The situation is similar for the convex hull formed ex¬ 
clusively by the pure product states. However since the 
product states is a manifold of measure zero in the set of 
all states, the convex hull of pure product states is located 
inside the convex hull of all pure states, apart from the 
points at the unit hypersphere corresponding to the pure 
product states. At the same time, the outer surface of 
the convex hull of product states does not separate posi¬ 
tive from negative matrices, and hence it must not exclu¬ 
sively contain degenerate matrices. In FigH] (a) we illus¬ 
trate the situation symbolically by showing pure product 
states as points at the spherical surface and the convex 
hull of these points by a polytope inside the sphere. At 
the surface and inside the polytope the states are sepa¬ 
rable. 

In FiglT](c) we illustrate that inseparable states are the 
mixed states inside the body symbolizing the convex hull 
of pure states but are outside the polytope symbolizing 
the convex hull of the product states. In FiglU (e)-(f) 
we illustrate the geometric meaning of Eq.®, that each 
mixed state can be represented as a sum of separable 
state symbolized by the polytope within a scaled sphere 
of radius a = entangled state in the body 

within a scaled sphere of radius b = bi = 1 — a. 

In the situation where b is minimum (6 = B, Eq.Q), 


the corresponding extremum states are at the surfaces 
of the polytope and the body, respectively, as shown in 
Figs[T](e)-(f). Obviously for a given a state such a decom¬ 
position is unique and the extremum state on the outer 
surface Sch of the convex hull corresponds to essentially 
entangled component in Eq.®. 

Let us now turn to the properties of the essentially en¬ 
tangled component pent which, as it will be shown now, 
is a density matrix of rank cIe strictly less than the di¬ 
mension N of the Hilbert space of the entire system. The 
essentially entangled component belongs to the outer hy¬ 
persurface Sch of the convex hull of all states, but not 
every state on Sch is an essentially entangled compo¬ 
nent; only some of them which do not contain the sepa¬ 
rable part (see FiglT](c)). In addition, the eigenvectors 
of Pent, I'*/’/) of Pent with 1 = 1,..., (Ie, are necessarily K- 
entangled pure states in the sense that these cannot be 
written as direct product of K pure states corresponding 
respectively to the K subsystems. Henceforth, we call 
pure states which are direct products of K pure states of 
the K subsystem, K-product states. 

Consider now the Hilbert space He of dimension dE, 
which is associated with the eigenvectors |'0i) of pent- 
Each state |'0) belonging to the Hilbert space He is ap¬ 
parently a linear combination of the eigenvectors, 

dE 

= ■ ( 8 ) 
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The vector can be also seen as a result of the action 
of an element Ue of the unitary group SU{(Ie) associated 
with the Hilbert subspace He sX one of the eigenvectors, 

\^) = IJeHi)- (9) 

Now, let us consider the convex hull of the states [Vi) of 
the subspace, which naturally contains pent- The condi¬ 
tion that Pent does not have any separable components, 
IV'prod) (V'prodI) implies that the convex hull does not con¬ 
tain a product state IV’prod) (V'prodI which is possible only 
if the Hilbert space He does not contain |^prod)- We 
name a Hilbert subspace with such a property an essen¬ 
tially entangled subspace of dimension dE and in what 
follows, with the help of this necessary condition, we find 
an upper bound on dE- 

Theorem II.1 The maximum rank dE-max of an essen¬ 
tially entangled component pent for a system of dimension 
N comprised by K subsystems each of them of dimension 
Nk, isN 

Proof Let us assume that the essentially entangled com¬ 
ponent Pent is a density matrix of rank dE, and consider 
the subspace He which is spanned by its if-entangled 
eigenvectors {\ijj\) ,\tp 2 ) ■ IV'dE)}- Let us also consider 

the orthogonal compliment of the subspace He, H^ of 
dimension N — and arbitrary select a set of mu¬ 
tually orthogonal vectors {|xi), |x 2 ), • ■ • \xN-dE)} span¬ 
ning H^. 

The subspace He is not essentially entangled, if there 
is at least one product state |'0prod) which can be ex¬ 
pressed as in Eq.®, 

dE 

l^/’prod) = l^i) • (10) 

where A’s are complex numbers. Equation m implies 
that |'0prod) must be orthogonal to every element {|xi)}i 
with i = 1,N — dE, oi the chosen basis in H^, 

(V^prodl Xi——ds) ~ 0 • (H) 

The maximum number of such conditions equals to the 
number of parameters defining a product state, which for 
a if-product state amounts to — K. Therefore 

the maximum rank of an essentially if-entangled density 
matrix cannot be equal or exceed — K. 

The maximum rank is smaller when we speak not 
about the essentially if-entangled component, but about 
the essentially entangled component which does not con¬ 
tain, not only if-product, but any product state. For this 
case one has to identify the bi-partition of the system that 
yields the maximum number of parameters characteriz¬ 
ing the product state. 

In the Appendix we provide a more detailed proof of 
this theorem. 


If this theorem is applied to the case of 2 qubits in mixed 
state, dEmax = I is obtained meaning that that the es¬ 
sentially entangled component can only be a pure entan- 
gled state. This result is in agreement with the results in 
[l^ where the bipartite case is treated. We would like to 
note, that the example studied in Section fVl gives some 
preliminary evidence that pent is prone to stay very near 
to pure states (Tr [pent] ~ 1) even though dEmax —t N 
ior N >> 1. 


III. THE LINEAR PROGRAMMING 
ITERATION ALGORITHM THAT YIELDS THE 
ESSENTIALLY ENTANGLED COMPONENT OF 
A DENSITY MATRIX 

One can find the maximum separable and the essen¬ 
tially entangled components of an arbitrary density ma¬ 
trix straightforwardly with the help of the linear pro¬ 
gramming algorithm applied to the convex hull of general 
pure states and the “polytope” of pure separable quan¬ 
tum states. The main obstacle on this way is a high di¬ 
mensionality of the corresponding Liouville vector space, 
which makes intractable the direct approach within any 
approximation. In fact, even for the simplest multipar¬ 
tite system of three qubits, the dimensionality (fV^) of 
the density matrix space is 64, such that even for the 
rather low-accuracy approximation attributing just 10 
points per dimension, one encounters a polytope of al¬ 
ready 10®^ vertices. 

Here, we suggest a way to crucially decrease the num¬ 
ber of the vertices that enter as samples in the algorithm 
and, in consequence, the computational complexity of 
the procedure. We first notice the fact that the solution 
of the problem and, in general, any convex decomposi¬ 
tion of the form Eq.([3]), allows for at most non-zero 
coefficients and bi- This observation can be formally 
justified by a theorem of Caratheodory as mentioned in 
Q. In the limit B = states are the 

vertices associated with the corners of the facets corre¬ 
sponding to the solution, as illustrated in Fig. [T] (e)-(f), 
while other vertices can be discarded. 

Therefore, at first step we may randomly take prod¬ 
uct states, N'^ general states and in order to ensure the al¬ 
gorithmic stability, complement this set by the eigen¬ 
vectors of the given density matrix. We find the solu¬ 
tion of the linear programming problem, which typically 
has complexity ^ (2Y‘^) , and thereby identify at most 
Y^ — J product states and J general states with nonzero 
coefficients ai and bi, respectively. The linear constraint 
imposed on the algorithm is the minimization of 
and the solution provided is a ‘local’ minimum, for the 
given set of vectors fed to the algorithm. Our aim is to 
find the global minimum value of ^^at is equal 

to B and to this end we create an iterative optimization 
loop which guides us there. 

At the second and subsequent steps, we take the 
product states resulting from the solution of the op- 
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timization problem at the former step and by apply¬ 
ing to each of them randomly chosen local trans¬ 
formations exp {i we generate ~ N'^ 

new product states. We also generate new entan¬ 
gled states by applying random generic transformations 
exp |i to each of the entangled states ob¬ 

tained at the former step. Here i numerates generators of 
the SU(N) group while i G local mean generators of the 
subgroup of local transformation. Random parameters 
and are normally distributed with width gradually de¬ 
creasing with the number of the iteration step. We again 
solve the linear programming problem for ^ N'^ vertices 
at these two new polytopes and iteratively repeat all the 
procedure till the result converges. Note that each next 
step, the presence of the solution of the former step of 
the loop is essential in order to guarantee an outcome 
from the linear programming algorithm. The set of the 
eigenvectors of the density matrix plays this role for the 
first step. Numerical inspection shows that the final re¬ 
sults of the algorithm i.e., the product component psep 
and the essentially entangled part Pent, Eqs. ®-© , are 
always the same for different runs. 

The algorithm described above concerns the case of 
full separability of a state or else, the identification of 
the essentially RT-entangled component. The same steps, 
can be applied if we make a repartition of the initial sys¬ 
tem and consider L-separability of the state with L < K. 
Furthermore, if the set of separable states is enlarged to 
include other special classes of pure states e.g. states of 
the W class [l3 ; ttien one can apply the idea of the algo¬ 
rithm for revealing the classification of mixed multipar¬ 
tite entangled state as the one introduced in [l^ for three 
qubits. We would like to add here that for the specific 
case of three qubits in mixed state, a lot of progress has 
been recently made on the classification of entan glem ent 
via analytic criteria and efficient algorithms |l(^-|la|. 

Finally it is important to mention that linear program¬ 
ming scales polynomially with the dimension of the vec¬ 
tor space under consideration in the general case but not 
always -still a zero-measure of non-polynomial cases may 
exist. In consequence, the same additional ‘rule’ has to 
be applied to the proposed algorithm and the identifi¬ 
cation of the special cases where the algorithm becomes 
non-polynomial is an interesting open problem, not re¬ 
solved in this work. However, on a practical level even 
in this case, a small random variation of the initial den¬ 
sity matrix brings the problem back to a polynomially 
complexity. 


IV. SUGGESTIONS FOR CHARACTERIZING 
ENTANGLEMENT PROPERTIES OF THE 
ESSENTIALLY ENTANGLED COMPONENT 


We may claim that all information relevant to entan¬ 
glement is contained in the essentially entangled part pent 
of the density matrix. Though this is not the main ob¬ 
ject of this work, we make some simple suggestions for 


analyzing ent ang lement properties of pent employing pre¬ 
vious results about characterization of entanglement 
for pure states. 

For pure quantum states, entanglement is directly re¬ 
lated to the factorizability at state vectors, and therefore 
one can characterize entanglement by identifying the or¬ 
bit of local transformations for a given state. This orbit 
can be marked by a complete set of polynomial invariants 
or alternatively by the coefficients {/3} of the tanglemeter 
({/?}) = of ^ state [T^. 

The state dehned as |c({/3})) = exp Nl ({/3}) |0) is the 


so called canonical state, and this can be reached from 
the state under study by the action of local operations 
under the constraint that the population of the reference 
state |0) is maximized. In addition to the identification 
of the orbit of local transformations the tanglemeter gen¬ 
eralizes the concept of logarithm to vectors and its coef¬ 
ficients straightforwardly reveal the factorization proper¬ 
ties of the state. 

Entanglement of mixed states cannot rely only on one 
operation of group multiplication but it should also in¬ 
volve the procedure of casting in convex sums. There¬ 
fore the algebraic structure does not suggest a natural 
framework for the characterization of entanglement in 
this case. Construction of an approach to entanglement 
characterization is a convenience just complementing the 
exhaustive information contained in the essentially en¬ 
tangled part of the density matrix. 

A straightforward way to characterize entanglement of 
mixed states would be to find the tanglemeters of the 
eigenstates of Pent- However, it does not mean that an 
entangled state corresponding to another orbit cannot 
be detected. In fact, any pure state which belongs to the 
essentially entangled subspace He spanned by the eigen¬ 
vectors of Pent 7 is also a legitimate representative of the 
ensemble of entangled states associated with this density 
matrix. One therefore may want to hnd the “corners” 
of this ensemble of states by identifying the state |ci) 
in He closest to the set of product states V, followed 
by identification of a state |c 2 ) T |ci) closest to V then, 
jca) T |c 2 ), |ci) etc, till \cdE)j and calculate tanglemeters 
for these “corners”. Tanglemeter coefficients of any state 
in He will therefore be within the borders given by these 
“corners”. We would like to mention here that the use of 
the tanglemeter as a method for characterizing multipar¬ 
tite entanglement is not essential here. One may apply 
this idea to other measures of multipartite entanglement 
for pure states as are the entanglement monotones from 
anti-linear operators introduced in 

One more option is to hnd the tanglemeter coefficients 
distribution function 


Pm) = /(C({;9(x)})|pe„t|c({/3(X)})) 

5{{l3{x) - |3})d^J.,^(zHE ( 12 ) 

resulting from averaging over the Haar measure ^Ix^He 
in the subspace He in according with the Eq.® with 
weight suggested by Pent 7 the probability to have canonic 
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state with given tanglemeter coeficients. The number 
P ({/3}) gives the probability density to find an entan¬ 
gled state which belongs to the orbit characterized by 
the set {/?} of the tanglemeter coefficients. In the case 
where one of the eigenvalues of pent is much larger than 
others, the probability distribution P ({/?}) locates near 
the tanglemeter of the corresponding eigenvector and it 
can be adequately characterized by a small covariance 
matrix of the tanglemeter’s coefficients. 


V. EXAMPLE 


We now present an illustration of the introduced meth¬ 
ods at a physical example of an open system experiencing 
loss of coherence due to presence of classical noise. The 
model comprises three elements: two two-level systems 
and a three-level system. The local symmetry group 
for each of two-level systems is the SU{2) group, for 
the three-level SU{3) group, while for the total assem¬ 
bly the group of transformations (local and non-local) is 
the SU(12) . We consider the following physical ingredi¬ 
ents of the combined system: an atom in p-state (L = 1, 
Ml = 0, ±1) in a static magnetic field, which paramet¬ 
rically interacts with a two-mode electromagnetic field. 
We also assume that each of the field modes allows for 
two possible polarizations of the photons. 

The Hamiltonian of the system consists of four parts: 

(i) the Hamiltonian of the first field mode Hi = 
(ajua; + apiy) with wavevector and polarizations 

X and y, 

(ii) the Hamiltonian of the second mode H 2 = 
kx (ppy + bp 2 '^ with wavevector kx, 

(iii) the Hamiltonian of the atom H 3 = in the 

static magnetic H = {Hx, Hy, H^} field, where L the an¬ 
gular momentum vector operator, and 

(iv) the Hamiltonian describing the parametric interac¬ 
tion 

^ (at a, + at a.) Xf , + Wy) YZ 

Ha = --^-r-: (Id) 

kz - kx- UJ 2 


which results from the second order perturbation the¬ 
ory applied over the dipole interaction (at -|- a^) X + 

PI + ay)Y+(pl + hy)Y+ {bl + bz)z. 

Here a] and 6| are the photon creation operators of 
the first and the second mode, respectively, correspond¬ 
ing to polarization along the direction i, while at and 
bi are their conjugate photon annihilation operators, re¬ 
spectively. By uji and a ;2 we denote the frequencies of the 
allowed dipole atomic transition from the state p, that are 
closest to the frequencies of the first kz and the second kx 
photon modes, respectively. The atomic optical electron 
radius-vector operator R = |a,F, and the angular 

momentum vector operator L = < L^,, Ly, > enter the 


Hamiltonian as the tensor product and the scalar prod¬ 
ucts with the magnetic field, respectively, while the light 
velocity, the electron charge, and the Planck’s constant 
are set to unity. 

Since parametric interaction implies conservation of 
the total number of photons on the two modes, H 1 +H 2 is 
an integral of motion for the system and only the Hamil¬ 
tonians H 3 and H 4 are responsible for the dynamical pro¬ 
cess of interest. The relevant part H = H 3 + H 4 can be 
re-written in a more convenient way, noting that the x, 
y, and 2 : components of the vector-operator L form an 
su{2) subalgebra of the symmetry algebra sm(3) of the 
atomic triplet p, while the operators XY and YZ enter¬ 
ing H 4 as the tensor product of the components of R do 
not belong to this subalgebra and yield other generators 
of SU{3) group. All these operators can be expressed 
in terms of Gell-Mann matrices Xi with i = 1,...,8. 
Moreover, the properly normalized bi-lineal photon op¬ 
erators alpy + apix, — apx, and alox — apy of the 
first mode form an su{ 2 ) algebra, and so do the similar 
operators of the second mode. Therefore, these can be 
expressed in terms of the Pauli matrices ai^i and ^ 2 , 1 , 
respectively, with i = x,y,z. Summarizing, the Hamilto¬ 
nian H = H 3 + H 4 can be caste in the form 

3 

H = E Xifi + f40'i,xXA+fe0'2,xXe+eiai^z+£20'2,z, (14) 

i=l 


where the parameters /i=i, 2,3 depend on the static held, 
parameters /4 and fe are governed by the detuning of the 
photon frequencies from the atomic transition frequen¬ 
cies, and parameters £ 1 = 1,2 deviate from zero when the 
photon frequency turns to be dependent on the polariza¬ 
tion in the presence of an anisotropicity of the refraction 
index (that is when kz is slightly different for the x and 
y polarizations, and similar for kx)- 

Now let us consider a realistic situation where the 
static held experiences small and rapid huctuations, that 
is fi{t) = fi + bfi{t) for i = 1,2,3. In this case the Li- 


ouville equation ip = 



describing the time evo¬ 


lution of the density matrix p{t) of the assembly, can be 
averaged over these rapid huctuations Sfi{t), yielding [Hj 
the following Lindblad master equation 



— ^ 


smsfpt) 




>^j,P 


(15) 


where the upper bar denotes time average. Substitution 
to this master equation in the Liouville representation 

143 

p{t) = '^nit)9r (16) 

i=0 

of the density matrix in terms of the generators of the 
unitary group SU (12), yields a system of 143 linear, hrst- 
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order differential equations 


irk 


Tlk 




143 


E (Tr{l?f 

771=1 


‘^1 

^ (17) 

3 

* j=i 


! 9m 

]} 


for the real vector components Vi (t). Straightforward an¬ 
alytic solution of Eq. (ED gives oscillations with time for 
some of the coefficients ri{t) while others die off with 
rates determined by the relaxation operator TZk,m- 
A considerable amount of work on the understanding 
of the dynamics of entanglement has been performed so 
far and we refer an interested reader to Q for a complete 
review and reference list. In Fig. [5] we graphically repre¬ 
sent a generic solution for this example, as a spiral in the 
Liouville space, gradually approaching a stationary solu¬ 
tion. This picture also provides a complementary point 
of view on the phenomenon of sudden death and revival 
of entanglement . With the course of time, we expect 
the essentially entangled part to oscillate between differ¬ 
ent subspaces and eventually to vanish for sometime - 
when the density matrix is passing inside the polytope 
of separable states as it is illustrated in Fig. The re¬ 
vival of entanglement is marked by the exit of the density 
matrix from the polytope. This graphical representation 
is justified by the calculations which we present in the 
following. 


c = 0.03 
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FIG. 2: A symbolic description of the trajectory in the Liou¬ 
ville space of a mixed state undergoing loss of coherence due 
to interaction with the environment. Crossing of the polytope 
of the separable states results in sudden death (or birth) of 
entanglement. In the inlet we list the numerical values of the 
parameters of the model. 

We now solve the model Eq. ED for a set of given val¬ 
ues for fi presented in Fig. [51 and reconstruct the den¬ 


sity matrix p{t) with the help of Eg. (1161) . We summa¬ 
rize the results of calculations in Fig. [3] In Fig. [3] (a) 
we plot the purity P{t) = Tr [p^(t)] of the density ma¬ 
trix as a function of time. At each time step we apply 
the algorithm and we decompose the density matrix as 
p(t) = (1-B(<))psep(t)+B(t)pent(t), Eq.dlD. InFig.|3](b) 
we plot the weight B{t) = Eq.dSD, of the essen¬ 

tially entangled component in the density matrix. The 
weight B{t) is decreasing with time faster than the pu¬ 
rity does and in addition it exhibits some oscillatory be¬ 
havior that can be probably explained by motion of the 
essentially entangled component along the facets of the 
polytope. In Fig. [3] (c) we plot the rank ds of pent(t), 
and we observe that this moves in a rather random way 
between the value 1 and 5. We note that if full {K = 3) 
separability is considered then dE max = 7. However in 
our program we have included in the “polytope” of sep¬ 
arable states also the bi-separable states, thus actually 
dEmax = 5. The “jumps” of the rank demonstrate the 
recursive move of essentially entangled component be¬ 
tween different essentially entangled subspaces on Sch- 
Moreover, in the time interval [18.8—19.7], B{t) vanishes 
implying that the state enters inside the polytope of sep¬ 
arable states. This physical situation describes a sudden 
death and sudden revival of entanglement a phenomenon 
[2^-[23 which has been studied extensively with other 
methods. Our geometric decomposition offers additional 
information on the origin of this phenomenon, see Fig. [2] 
In order to analyze the entanglement properties of the 
essentially entangled component we first note that for 
the chosen model system in the vast majority of the time 
steps, there is a dominant eigenvector Cdom for Pent with 
a corresponding eigenvalue Xdom > 0.9, see Fig. |3| (d). 
Therefore, for this specihc example and assigned param¬ 
eters, it makes sense just to analyze entanglement prop¬ 
erties of edom 7 whenever the condition Xdom > 0.9 is satis¬ 
fied, and to conclude from this analysis the entanglement 
properties of Pent- Naturally, this analysis together with 
the weight B{t), give all the information necessary to 
describe entanglement in p. 


We analyze the entanglement properties of edom with 
the help of the method of nilpotent polynomials M- In 
the Appendix we provide an explicit method for deriving 
the general expression for the tanglemeter of a wavevec- 
tor describing an assembly of a three-level system and 
two two-level systems: 

a2 +/doiidt ^2 + 
l32iou'^(rt + I320 iu~''^2 + 

I3iiii+atd+) (18) 

with /3iii,/3201)/3210)/3iio being positive numbers and 
PioijPoii being complex. The matrix representation of 
the nilpotent variables (operators) u~^, v~^, (T+ is also 
provided in the Appendix. Concerning now the physical 













FIG. 3: We solve the Lindbland equation for the example and we apply the algorithm at each time step, (a) Purity of the 
assembly, (b) The statistical contribition of Pent(t) to the density matrix, (c) The rank of Pent{t). (d) The eigenvalue of 
the dominant eigenvector of Pent(t). (e)-(h) The oscillations of the real coefficients of the tanglemeter. In the time interval 
[18.8 — 19.7] sudden death of entanglement takes place and then its revival. 


meaning of the coefhcients. The coefficients of the tan¬ 
glemeter even though are not entanglement monotones 
3 in the strict sense, these are invariant under the action 
of local transformations and the presence of any non-zero 
term in the tanglemeter ensures the presence of entan¬ 
glement. More precisely, the coefficient fim ensures the 


presence of genuine tripartite entanglement in the state 
while the rest of the coefficients are related to bipartite 
entanglement. In Fig. [3] (e)-(h) we plot those coefficients 
which are positive, and we observe that these oscillate 
with time without dissipation. The same holds for the 
real and imaginary parts of the complex ones not shown 
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on the figure. 

With this example, in addition to the death and re¬ 
vival of entanglement, we observe two interesting phe¬ 
nomena which need more case study in order to decide 
whether are specific to this example or general. The first 
is the presence of a dominant eigenvector in the essential 
entangled component and the second is the oscillations 
without dissipation of the entanglement characteristics 
of the essential entangled component. 


VI. CONCLUSIONS 

In this work we have studied a concept related to en¬ 
tanglement of mixed states namely the essentially entan¬ 
gled component of a mixed multipartite state and more 
important, we have suggested an efficient algorithm for 
its identification. The essentially entangled component is 
the complementary part to the best separable approxima¬ 
tion introduced in [l2| and this naturally contains all the 
entanglement of the density matrix. We analyze some 
properties of the essentially entangled components and 
we suggest methods for characterizing its entanglement 
content. 

Our main tool is the accustomed geometric description 
of mixed quantum states in the spirit of Bloch vector rep¬ 
resentation, which results from the decomposition of a 
density matrix over the generators of the relevant group. 
We have shown that pure states are not everywhere on 
the surface of this hypersphere, in the contract to the 
Bloch vector, and that the convex hull of pure states 
from a convex “body” inside the sphere. The convex hull 
of separable states forms a convex “polytope” inside the 
“body” of general states. As a consequence the entan¬ 
gled states inside the body and outside the polytope can 
be represented as sum of a separable state on the surface 
of the polytope and an essentially entangled component 
located on the surface of the “body”. This geometric 
picture gives the guidance for constructing the algorithm 
and for analyzing the properties of the essentially entan¬ 
gled component. The latter being located on the surface 
of the “body”, form there sets of lower dimensions, such 
that the rank of the relevant density matrix does not ex¬ 
ceed a number which depends on the dimensions of the 
total system, and on its chosen partition. 

Finally, at a particular example we study the dynam¬ 
ics of an open quantum system and we reconstruct the 
time trajectory of the decomposed density matrix inside 
the convex “body”. Sudden death and sudden birth of 
entanglement can be seen as the results of crossing of the 
of the trajectory of the density matrix with the surface 
of the “polytope” of separable states. There are some 
other interesting phenomena appearing in this example 
but these still need further studies to lead to general con¬ 
clusions. 

Concerning possible applications of the results. The al¬ 
gorithm introduced in this work scales polynomially with 
the dimension of the system in the general case, and it 


can be employed to study open questions about entangle¬ 
ment in mixed states. For instance, this can be applied 
straightforwardly to address the question of the relative 
volume of separable states over entangled mixed states as 
function of the total purity of the system and the total 
dimension of the system [25l | . An answer to this question 
can serve to the evaluation of emerging quantum tech¬ 
nologies and their quantum limits. Moreover, the essen¬ 
tially entangled component containing all entanglement 
properties of the density matrix may also provide new di¬ 
rections to entanglement detection |26l | and entanglement 
distillation techniques. 
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APPENDIX 

A. A second formulation and proof of the main 
theorem 

Here we provide a more detailed formulation and proof 
for the theorem given in Section HIl which does not relay 
on a particular quantum mechanical representation. 

The maximum rank dEmax of an essentially entangled 
component is Nca — Ncs, where Ncg is the dimension 
of the Cartan subgroup of the group of all transforma¬ 
tions on the state and Ncs is the dimension of Cartan 
(sub)subgroup generating only local transformations. 

Remark. The numbers Ncg and Ncs give the num¬ 
bers of complex parameters characterizing generic and 
product state vectors, respectively, on V = Ncg + 1 di¬ 
mensional Hilbert spaces. 

Proof Consider a density matrix p and its decompo¬ 
sition to the essentially entangled and separable part 
p = (1 — B)psep + Bpent- Since B corresponds to a min¬ 
imum value of all possible weights, we conclude that no 
e > 0 and product vector \p) exist such that Pent —e \p) {p\ 
is a positive matrix. Considering now the essentially en¬ 
tangled subspace He spanned by the eigenvectors |'0i) 
with non-zero eigenvalues of pent with i = 1 ,... ,dE, this 
condition means that no product state \p) exists in He. 
Indeed, for the case where 

dE 

\p) = ^\f^^) + e'\p') (A-1) 
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with (p'l ij^i) = 0 for every i = 1 ,..., dg one identifies the 
vector Ip') orthogonal to the subspace of eigenvectors 
which makes 

(p'Kdent - eb) {p\)\p) 

= -e|(p|p')|"<0, (A-2) 

and therefore extremality implies that no product state 
is orthogonal to the orthogonal compliment ot He 
spanned by the eigenvectors \ipi) of Pent with zero eigen¬ 
values and i = cIe, • ■ ■, Ncg- 

In other words, in order to find such a state we have 
to satisfy Ncg — dE equations (p| ipi) = 0 with f = d_E-|- 
1 ,..., Ncg for a product state |p) given by specification 
of its Ncs parameters. This is impossible when Ncg — 
dE ^ Ngs, which determines the maximum rank dE max 

of Pent- 

B. Deriving the tanglemeter of the physical 
example in Section [Vl 

The system under consideration consists of the two 
modes of the held interacting with a three-level atom. 
The Hilbert space thus is of dimension iV = 12, a direct 
product of the spaces of two two-level systems (qubits) 
and of one three level system (qutrit). In the standard 
computational basis a state vector of the system is ex¬ 
pressed as 

1^) = boDO 1000) -I- bioo 1100) -I- '0200 1200) -|- 0oio 1010) -I- 
00011001) 0110 |110) -I- 01011101) -I- 00111011) -I- 
0210 |210) -I- 0201 |201) -I- 0111 |111) -I- 0211 |211). 

or alternatively using the nilpotent creation operators 

/O 0 1\ 

u+ = 0 0 0 , (A-3) 

\o 0 oy 


/o 0 o\ 

t+ = 0 0 1 , (A-4) 

\o 0 oy 



as 

1^) = (0000 + 01OO1’'' + i'200U+ + 0oioo-y + 00010'^ 

i!iwi+cF^ + iiimi+a^ +^2wu+d^ + 

02 OlU'''d^ + ^2 + 021 lU+(Ty(J^) | 000 ). 

The next step that should be performed is the applica¬ 
tion of all the available local transformations {SU 
1, l(8)S'C/(2)® 1, 1®1®S'{7(2)) on the given state Idi) in 
order to construct the corresponding canonic state |d>c) 
which marks the orbit of local transformations. To sim¬ 
plify the procedure, we apply the local transformations 
on a given lih) in the following order: 

(a) We hrst apply local operations generated by the 
operators |(Ti , df, df, df, A 4 , A 5 , Ae, Arj and we require 

that the polulation of the reference level | 000 ) is getting 
maximum. Under this condition the populations of the 
levels : | 100 ), | 200 ), | 010 ), | 001 ) are vanishing. 

(b) We then apply local operations generated by 
Ai,A 2 | to maximize also the population of the level 

I 111 ). This way the contribution of the level |211) also 
vanishes. 

(c) Finally we apply local operations generated 
by |(Tf, df, A 3 , Asj in order to make the phase of 

I 111 ), | 210 ), | 201 ), 1 110 ) equal to the phase of the am¬ 
plitude of the reference level | 000 ). 

After this procedure one obtains the following form for 
the unormalized canonic state: 

l^'c) = (1 + aiioi+d^ + aioit+dy 

+aoiid^^2 o:2iou^ -I-a 2 oid+d^ 

-baiiit+d+d^) |000) (A- 6 ) 

with oiii, a 2 oi, 0 : 210 : Q^iio being positive numbers and 
0101,0011 being complex. 

The final step for arriving to the tanglemeter Nl {{^}) 
of the state is to take the logarithm of the polynomial on 
the nilpotent variables £+, d ^2 Ea. (IA- 6 l) . It is easy to 

show that 

({/ 3 }) = Piioi+df + (3ioii+d^ + Poud^d^ 
02iod+dy +P20iu^d} + I3iiii+d^d^ 

with 0110 = Olio, 0101 = oioi etc. 
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